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ABSTRACT 


High range resolution radar systems have many advantages such as target 
classification, resolution of multiple targets, accurate range measurement, target range 
profile and detection of low radar cross section (RCS) targets in clutter. High range 
resolution requires large bandwidths. Stepped frequency waveforms can achieve high 
range resolution by increasing the effective bandwidth without increasing the 
instantaneous bandwidth which would increase the hardware requirements including 
higher analog to digital (A/D) sampling rates which are limited by existing technology. 
Under today’s hardware limitations, the stepped frequency waveform becomes very 
important. This thesis briefly discusses the stepped frequency radar and associated signal 
processing, it investigates the ambiguity function of the stepped frequency waveform and 
the stepped frequency radar system. Mathematical expressions of ambiguity functions are 
derived and the improvement of clutter suppression capability for the stepped frequency 
radar by rejecting initial pulses is also discussed. 
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I. INTRODUCTION 


High range resolution radar systems have many advantages such as target 
classification, resolution of multiple targets, accurate range measurement, target range 
profile and detection of low radar cross section (RCS) targets in clutter. High range 
resolution is typically achieved either by decreasing the duration of the transmitted signal 
or by modulating a transmitted signal of relatively longer duration. In either case, the 
instantaneous bandwidth of the signal goes up and thereby increases the hardware 
requirements to include higher A/D sampling rates (which are limited by existing 
technology). However, with the stepped frequency waveform, it is possible to achieve 
higher range resolutions with lower instantaneous bandwidths and lower A/D sampling 
rates. Irrespective of the waveform or compression method used, high range resolution 
requires large bandwidths. For the stepped frequency waveform, a large bandwidth is 
obtained sequentially by changing the carrier frequency in steps over several pulses instead 
of within a single pulse. This waveform potentially can be implemented in high 
performance radar systems with lower A/D sampling rates, which is one of the current 
bottlenecks in the development of radar systems. Because of the importance of the 
stepped frequency waveform, this thesis briefly discusses the stepped frequency radar, its 
associated signal processing, and analyzes it by the ambiguity function. In Chapter II the 
stepped frequency waveform and its corresponding radar system are introduced. In 
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Chapter in the matched filter and the correlation definition of the ambiguity function are 
first discussed. The properties and applications of the ambiguity function are also briefly 
discussed. The mathematical expression of the auto-ambiguity function for the stepped 
frequency waveform is then derived from the definition and is verified by comparing the 
ambiguity diagram with one obtained by simulation. Finally, the mathematical expression 
of the cross-ambiguity function for the stepped frequency radar signal processor is derived 
and verified, and the improvement of the clutter suppression capability for the stepped 
frequency radar by the rejection of initial pulses is discussed. In Appendix A, the 
ambiguity functions for the single pulse, constant frequency pulse train, linear frequency 
modulated (LFM) pulse, and the discrete frequency modulated pulse are given for the 
purpose of comparison. Mathematical expressions of the ambiguity functions for these 
waveforms are derived and verified. In Appendix B, the computer programs of ambiguity 
functions for all the waveforms discussed in the thesis are listed. 
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H. OPERATING PRINCIPLES OF THE STEPPED FREQUENCY RADAR 

A. STEPPED FREQUENCY WAVEFORM 

The stepped frequency waveform can be described by a series of N coherent 
pulses, where the carrier frequency increases from pulse to pulse by a fixed increment A f 
as shown in Figure 1. Each pulse has a fixed pulse width T s , and the pulses are transmitted 
at a fixed pulse repetition interval (PRI). The frequency of the kth pulse is given by 

f k =T 0 +(k-l)A/, (1) 

where f 0 is the nominal carrier frequency and A/is the frequency step size. [Ref. 1, pp. 
160-161] 

The instantaneous bandwidth of this waveform is approximately equal to the 
inverse of the pulse width and is much less than the effective bandwidth. The waveform's 
effective bandwidth, denoted as B eff , is determined by the product of the number of pulses 
N, and the frequency step size A/as given below. 

Brff=NAf (2) 

The range resolution for any waveform is dependent on the effective bandwidth of the 
waveform. The range resolution of the stepped frequency waveform is given as 

Ar= 2^ (3) 

where c is the speed of the light. Thus the range resolution can be made finer by increasing 

the effective bandwidth NAf. [Ref. 2, p. 234] 
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1 2 3 4 N 

fO fO+Af fD+2Af fD+3Af fO+(N-l)Af 



N : number of pulses Tr : pulse repetition interval (PRI) 

Ts : pulse width Af; frequency step size 

fO : nominal carrier frequency 

Figure 1. The Stepped Frequency Waveform 


4 




The effective bandwidth is achieved by the signal processing of the data obtained 
during the coherent processing interval (CPI), which usually is the time on target. The CPI 

is given by the product of the number of pulses N and the pulse repetition interval (PRI). 

CPI=NxPRI (4) 

B. SYSTEM DESCRIPTION 

The block diagram of a radar which employs a stepped frequency waveform is 
shown in Figure 2. The transmitted signal in this radar is generated as in conventional 
coherent radars by mixing RF frequency (STALO) and IF frequency (COHO). However, 
the frequency step in each pulse is added by the frequency synthesizer as shown in Figure 
2. On the receiver side, a reverse operation of taking out the RF, step, and IF frequencies 
is performed. The IF amplifier serves as a matched filter to the envelope of each individual 
pulse. The output of the synchronous detector is a video signal in the form of in-phase and 
quadrature (I-Q) components. The output of the A/D is stored and organized according to 
range bins. Received signals from all N pulses in each range bin is collected. This data may 
be weighted before taking the DFT. The output of the DFT represents a finely resolved 
range bin. The size of the original range bin is cx/2. However, after DFT, the width of the 
subdivision of a finely resolved range bin is cx/2NAfx, which is equivalent to having a 
bandwidth of NAf [Ref.2, p. 236]. The unambiguous range R^ for this waveform is c/2Af. 
This will be quite short but it may not be a problem as it is applicable within the same 
range bin. For avoiding wrap around the target, R u should be greater than the target 
extent. 
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Figure 2. System block diagram for the Stepped Frequency Radar 
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m. ANALYSIS OF THE STEPPED FREQUENCY RADAR 
BY THE AMBIGUITY FUNCTION 

A. INTRODUCTION TO AMBIGUITY FUNCTION 

1. Matched Filter 

The matched filter maximizes the output peak signal to noise power ratio for 
the detection of a known signal. In the presence of additive white Gaussian noise, the 
impulse response of the matched filter is the complex conjugate of the image of the 
received waveform, h(t)=s*(t r t), where s(t) is the received signal, that is, it is the same as 
the received signal run backward in time starting from fixed time t r [Ref. 3, p. 261] In 
order to simplify the mathematical developement, it is customary to set tj at zero. The 
output of the matched filter is then written as 

oo oo 

y(t) = f s(k)h(t- X)dk = J s(X)s*(k - t)dk (5) 

—-oo —oo 

2. Definition of the Ambiguity Function 

If the transmitted signal S t (t) is u(t)e j27 ^ ot , then received signal S r (t) is given 

by u(t~To)e J ' 2n ^ 0+ ^ < ' t ~ TD) , where To is the unknown round trip time delay given by ^r- 
and fd is the Doppler shift of the returned signal. Assuming that the radar receiver is 
matched to the transmitted signal, we obtain the output of the matched filter as 

So[t d ) =2 u(t - To xu*(t-T D yM‘- T >)dt , ( 6 ) 

/ 

where T D is the estimate of time delay. 
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/ / 

It is customary to set T D and fo equal to zero, and to define Td~T d as x. The 

output of the matched filter is then 

So(t) = e~ j2 ^ dZ x J u(t-x)u*(t)e j2 ^ dt dt. (7) 

—oo 

Since the interest is only in the magnitude of the matched filter output, 
neglecting the first exponential term, we obtain the output of the matched filter as 

X(x,f d ) =J u(t-z)u*(t)e^‘dt (8) 

—OO 

In this form a positive x indicates a target beyond the reference delay t 0 , 
and a positive / d indicates an incoming target [Ref. 4]. In the literature the term 
"Ambiguity Function" is interchangeably used for X(x/ d ), IX(x/ d )l, and IX(x/ d )l 2 . In this 
thesis, the plot of IX(x/ d )l 2 is called the ambiguity diagram. If the receiver filter is matched 
to the transmitted signal, IX(x/ d )l 2 is termed the auto-ambiguity function, otherwise, it is 
termed the cross-ambiguity function. 

3. Properties of the Ambiguity Function 

The function \X(x,fd)\ 2 has the following properties [Ref. 5, p. 412]: 


Maximum value of \X(x,fd)\ 2 = lx(0,0)l 2 = (2 E) 2 , 

(9) 

\X(-x,-f d )\ 2 = \X(x,f d )\ 2 , 

(10) 

IX(t,0)I 2 = I J u{t-x)u*{t)dt\ 2 , 

—OO 

(11) 

\X(()J'd)\ 2 = 1 J u 2 (t)e j27 ^ d dt \ 2 , 

—oo 

(12) 
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U \X(xJ d )\ 2 dxdf d = (2 E) 2 . (13) 

The first equation given above. Equation 9, states that the maximum value of 

the ambiguity function occurs at the origin and its value is (2E) 2 , where E is the energy 
contained in the echo signal. Equation 10 shows that the ambiguity function is symmetric 
about the origin. Equation 11 describes that the ambiguity function along the time delay 
axis is the auto-correlation function of the complex envelope of the transmitted signal. 
Equation 12 shows that along the frequency shift axis the ambiguity function is 
proportional to the spectrum of u 2 (t). Equation 13 states that the total volume under 
the ambiguity function is a constant equal to (2E) 2 . 

The ideal ambiguity diagram consists of a single peak of infinitesimal thickness 
at the origin and is zero everywhere else, as shown in Figure 3. The single spike eliminates 
any ambiguities, and its infinitesimal thickness at the origin permits the frequency and the 
echo delay time to be determined simultaneously to as high a degree of accuracy as 
desired. It also permits the resolution of two targets no matter how close together they are 
on the ambiguity diagram. Because of two restrictions the ideal ambiguity function is not 
possible; first, the maximum height of the ambiguity function must be equal to (2E) 2 , 
which is limited, and second, the total area under the ambiguity function must be finite 
and equal to (2E) 2 . A reasonable approximation to the ideal ambiguity function might 
appear as shown in Figure 4. This waveform does not result in ambiguity since there is 
only one peak, but the single peak might be too broad to satisfy the requirements of 
accuracy and resolution. If the peak is made too narrow, the requirement for a constant 
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Figure 3. Ideal ambiguity diagram "after Ref. [5], 
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volume might cause peaks to form at regions of the ambiguity diagram other than the 
origin and give rise to ambiguities. Thus the requirements for accuracy and ambiguity may 
not always be simultaneously possible to satisfy. [Ref. 5, pp. 412-413] 

The particular waveform transmitted by a radar is chosen to satisfy the 
requirements for (1) detection, (2) measurement accuracy, (3) resolution, (4) ambiguity, 
and (5) clutter rejection. The ambiguity diagram may be used to assess qualitatively how 
well a waveform can achieve these requirements. The maximum value of the ambiguity 
function is an indication of the radar’s detection capabilities. Since the ambiguity diagram 
is often normalized so that the maximum value is equal to one, the ambiguity diagram is 
seldom used to assess the detection capabilities of the waveform. The accuracy with which 
the range and the velocity can be measured by a particular waveform depends on the 
width of the spike, centered at IX(0,0)l 2 , along the time and frequency axes. The 
resolution is also related to the width of the central spike, but in order to resolve two 
closely located targets the central spike must be isolated. It cannot have any high peaks 
nearby that might mask another target close to the desired target. A waveform that yields 
good resolution will also yield good accuracy, but the reverse is not always true. The 
presence of additional spikes can lead to ambiguity in the measurement of target 
parameters. An ambiguity measurement is one in which there is more than one choice 
available for the correct value of a parameter, but only one choice is appropriate. Thus the 
correct value is uncertain. The ambiguity diagram may be used to determine the ability of a 
waveform to reject clutter. If the transmitted waveform is to have good clutter-rejection 
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properties, the ambiguity diagram should have little or no response in the regions of clutter 
[Ref. 5, pp. 418-420]. Some examples of the ambiguity function are shown in Appendix A 
which includes single pulse, pulse train, linear frequency modulation pulse, and stepped 
frequency modulation pulse. 


B. DERIVATION OF THE AUTO-AMBIGUITY FUNCTION FOR THE 

STEPPED FREQUENCY RADAR 


In order to derive the auto-ambiguity function, according to the definition, 

Equation 8, we need the expression of the complex envelope of the transmitted signal. The 

transmitted signal of the stepped frequency radar consisting of a train of N pulses and 

changing carrier frequency from pulse to pulse by fixed increment A/ is shown in Figure 
1. It is mathematically represented as 


S t (t) = a t X 1 u(t - nT r )e^o+nm ? ( 14 ) 

n=0 

where u(t-nT r ) = l, if nT r <t< nT r + T s , 

0, elsewhere, (15) 

fstep - nAf , n = 0, 1 ,and 

fo = fcoho ^ f state , 

fcoho : the frequency of the stable coherent oscillator, 
f state : the frequency of the stable local oscillator, 
fstep : step frequency, 

A/ : the frequency step size, 

N : the number of the pulse sequences, 

T r : pulse repetition interval (PRI), 

T s : pulse width, 

a t : the amplitude of the transmitted signal. 
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The transmitted signal can be rewritten as 


St(t) — at 


X 1 u(t - nT r )e J ' lm ^ t 

. n=0 


*Wo* 


= a t U{t)e^ t (16) 

AM 

where U(t) = X u(t - nT r )e Jl7lnAfl is the complex envelope of the pulse sequence. The 

n -0 

signal received due to a discrete scatterer can be written as 

S r (t) = a r U{t-x)eP si trtJtri (17) 

where 


a r : the amplitude constant, 

x : the round trip time between the radar and the scatterer, 

fd : the Doppler frequency shift due to the normal relative velocity between the 

radar and the scatterer, 

x = f, (18) 

and 

U = 2j rh- (19) 

Substituting the complex envelope, U(t) , of the transmitted signal into 

Equation 8, we obtain the auto-ambiguity function of the stepped frequency radar as 


X(x ,f d ) = J X 1 u(t—mT r - t)^h) x Y u t (t-riT r )e~ ilml ^ x e^dt. (20) 

—o©m=0 n=0 


/ 

Changing variable t-nT r by t , we change the above equation to 


N-l/V—i 

X(xjd) = X X e~ j2nn2Afrr x e j 2nmn ^ T r X e~ j2rLmAft x e j2 ^ dnTr 

m=0/t=0 


x J u* (t luf/- (m - n) - x dt . 


( 21 ) 
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Defining the round trip time T in terms of integer and fraction parts 


i = (p+y)T r , (22) 

where p is an integer, y is the fraction, and 0 < y < 1 , we obtain the ambiguity function as 

X(x,f d ) = (Ts - iTr) X e^-P^Ts^Tr) x e j2npAfi x 

nfj-py)(Ts-*tTr) 


min (V-l fl-l+p) 

: 2 gfiKnifj-pAfiTr x g—j2mAfz 

n= max (Op) 


for 0 <y<£, 


(23-1) 


: (-7 + c Y _ nr ^ x £MHp+i)4/)(7>(y-i)7V) x sinrc(/Wp+i)4m?>CH)7V) 
v VI ^ jc(^-(p+i)y)(r s +(r-i)r r ) 


min(iV-l jV+z?) 

y/ e j2n(p+\)AfZys £ ^/27w(/^-(p+l)A/)rr ^ g-jlnnAfi; 


for 1 <y< 1 , and 


^O, 


Tr 


elsewhere. 


(23-2) 


(23-3) 


After simplifying the summation term and taking the absolute value, we obtain the 
expression for the auto-ambiguity function as 

I v/v s \l -/'r „/r \ w I smK(f<r-p6f)(Ts-iTr) | w | smn(f d -(2p^/)jy)(N-\p\)T r | 

\X(Xj d )\ - (i s ~iT r ) Xx I. 1 ’ 


for 0 <y<^f, 


(24-1) 
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= 0, elsewhere. (26-3) 

One can easily show that the auto-ambiguity function of the stepped frequency 

radar will become the ambiguity function of the pulse radar (shown in Appendix A) by 

setting A/ to zero. The ambiguity diagram, that is the plot of IX(x/ d )l 2 for the stepped 

frequency waveform for the following parameters: number of pulses N equals to 4, 
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duty-cycle equals 0.2, and frequency step size Af equals 2PRF is plotted in Figure 5. This 
plot somewhat looks like the ambiguity function of the pulse train rotated at an angle. The 
contour plot of the ambiguity function is shown in Figure 6. Plots of a cut along the 
frequency axis for t= 0 and a cut along the time delay axis for/ d =0 are shown in Figures 7 
and 8. It can be seen from Figure 7 that the spikes appear at multiples of PRF in the 
frequency axis. The first null for each spike is l/(NT r ) away from the center of the 
spike. Therefore, each spike has null to null width of 2/(NT r ) along the frequency axis. 
Amplitudes of those spikes vary as a SINC function. Zero amplitude will appear at l/T s . 
From Figure 8, there is 2N-1 spikes along the time delay axis. From Eq. 26-1, it is easy 
to show that the first null for the central spike (p=0) occurrs at l/(NAf). Accordingly, if 
1 /(NAf) is less than pulse width T s , the stepped frequency radar will have better range 
resolution than the pulse radar. Plots of cut along the delay axis with different frequency 
step Af,s are presented in Figure 9. The larger the radar's frequency step, the better the 


resolution of the radar. 


o o 



frequency (PRF) ' 10 ' 4 delay (PRI) 


Figure 5. Auto-ambiguity Diagram of the Stepped Frequency Radar, 
Pulse Number N=4, Duty-cycle=0.2, A/=2PRF. 


17 















frequency (PRF) 



delay (PRI) 


Figure 6. Contour Plot of the Auto-ambiguity Diagram for the Stepped 

Frequency Radar, Pulse Number N=4, Duty-cycle=0.2 ,A/=2PRF. 
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- 5 0 5 

frequency (prf) 

Figure 7. Plot of Cut along the Frequency Axis for x = 0 of the Auto-ambiguity Function 
for the Stepped Frequency Radar, Pulse Number N=4, Duty-cycle=0.2 , Af= 



Figure 8. Plot of Cut along the Time Axis for f=0 of the Auto-ambiguity Function for the 
Stepped Frequency Radar, Pulse Number N=4, Duty-cycle=0.2, A/=2PRF. 
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C. DERIVATION OF THE CROSS-AMBIGUITY FUNCTION FOR THE 


STEPPED FREQUENCY RADAR 

In this section the mathematical expression of the ambiguity function for the stepped 
frequency radar receiver and processor as shown in Figure 2 are derived. Since the filter 
used in this system is only matched to the individual pulse instead of matched to the 
transmitted signal which is composed of N pulses, the ambiguity function of this stepped 
frequency radar processor is called the cross-ambiguity function. The stepped frequency 
radar can be mathematically modeled as shown ih Figure 10. The transmitted and received 
signals were defined in Equation 16 and Equation 17 respectively. The received signal is 
mixed with a local oscillator, stepped frequency synthesizer, and a coherent oscillator. The 
stepped frequency synthesizer is synchronized so that the transmitter and receiver are on 
the same frequency step. As a result, multiple time around echoes will have frequencies 
which vary by multiples of frequency step size A/. The signal is then passed through a filter 
which is matched to a single pulse envelope. The output of the pulse-envelope matched 
filter is sampled and the samples undergo the discrete-time signal processing shown in 
Figure 10. 

The output of the stepped frequency synthesizer used in the receiver is defined as 

Ss'epit) =i |i(i- nT r )e-J 2 ^ 1 , (27) 

«=0 

where 

\iit-nTr) = 1, if nT r <t<(n+ 1)77, and 

0 , elsewhere. (28) 
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The output of the matched filter, z(t), can be derived by taking the convolution of 
the impulse response of the matched filter, h(t), and the received signal at the input of the 
matched filter. The received signal at the input of the matched filter can be defined as 

Si(t) = e~ j2 ^ ot x Sstep ( t ) x S r (t ). (29) 



The above equation can be rewritten as 


z(to) =a r j e j2v ?° x x e ,2 ^ d(y ^ l) xJxh*(to~ T|)G?r|, 


where 


J = L (X(ti- m r r )e-^xX w(r]-x- nT r )eJ 2m ^^ 

m=0 rv =0 


(31) 

(32) 


and x is defined in Equation 22. 


If each pulse of the received signal is assumed to fall within the same pulse of 


the stepped frequency signal no matter how x varies, then J can be expressed as 


min (N-l-pfi-~l) 

J = e -j2KKm x u(T\-x-nT r )e- j2mA ^, 

n ~max (0-p) 


where k =p, if pT r <x< (p+ l)T r -T s , 

=p + 1, if (p + l)7Y-r,<T<(p+l)7Y. 


Substituting Equation 33 into Equation 31, we obtain the output of the matched 


filter as 


z(to) = dre-Wf' Wa)t x J eWd^Utox 


min (N-l-pfl—1) 

X u(T) -x- nT r )e- j2ltn ^ x h*(t 0 - T])dn. 

(0—p) 


Letting E, = q - to , we rewrite the output of the matched filter as 


z(to) = a r e~ j2n ^ +fd)x xe ' 2ltM(|1 


00 P mm(N-l-pJV-l) 

xj e }2l W<r-m£, x ufe + to-x-nTr) xh*(-Qd^ 
-w L (0 ,—p) J 


= a r e x e j2 ^ mt °x J e j2 ^ x V(E, -(x-to),x)xh* (-fydt,, 

—CO 

where f m =fd~KAf, 


min (AM -p AM) 

V(t,x)= X u(t-nTr)e~J 2m &. 

/i=max(0 p) 
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It is customary to set t 0 to 0 for the computation of the ambiguity function and 
neglect the first exponential term, since we are only interested in the magnitude of the 
matched filter outputs. Redefining the output of the matched filter z(t 0 ) as X(x,f m ), we 

obtain the output of the matched filter as 

X(z,f m ) = J x Vit - x, x) x h* ( -t)dt . (39) 

-OO 

According to the definition of the ambiguity function in Equation 8, we treat f m as 
the frequency shift of signal at the input of the matched filter. This frequency shift consists 
of the Doppler shift and the frequency changes due to the stepped frequency waveform. 
V(t,x) is then the envelope of signal at the input of the matched filter. Thus X(x/ m ) 
becomes the cross-ambiguity function of the stepped frequency radar without the 
processor. In this form a positive x indicates a target beyond the reference delay t 0 . 

In the rest of the section, the ambiguity function is extended to include the signal 
processor. The matched filter output is sampled at instants t=nT s and the samples 
undergo discrete-time signal processing. The sample rate, l/T s , is greater than or 
equal to the bandwidth of the pulse envelope u(t). Further, the pulse repetition interval 
(PRI), T r , is an integer multiple, k, of T s . The sample stream for all the pulses is sorted 
into k parallel processing streams termed ambiguous range bins. Processing within each 
range bin utilizes only samples from time instants separated by integer multiples of T r . 
[Ref. 6] 
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The ith sample for range bin k can be defined as 

z k (i) = z(kT s + iT r ). (40) 

These samples from the matched filter output which are associated with a range 

bin are processed together by a discrete Fourier transformation (DFT) processor. A 

Hamming window is applied before executing the DFT to reduce frequency sidelobes. 

Out of N samples in a range bin, the first I samples are rejected to allow for the settling of 

transients due to range ambiguous clutter [Ref. 7]. DFT is taken as the last M weighted 

samples of each range bin as given below 

Z(k, l) = L w(i)e~ j2ni Mz k (i +/), (41) 

i =0 

where 

k : range bin index, 
l : Doppler filter index, 

M : the number of samples that DFT processed, and M <N-\p\ , 

i : sample index, i=0,l,...,M-l, 
w : the window coefficients, 

I : the number of samples that will be rejected initially. 


Substituting Equations 35 and 40 into Equation 41, we obtain the /th output of the 

DFT for range bin k as 

Z(k,l) = a r e-J 2K ^ x x 


00 min (AM -p /VM) 

J 2 u(x\-x-nT^e- 11 ™^ x h*(kT s -r\)dT], 

—=o n -max (0~p) 

where 


M-l n . / 

K it) = I ( t+ (i +r)T r ). 


(42) 


(43) 
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/z’,(t) is the impulse response of the radar receiver processor which includes the 
pulse envelope matched filter, weighting function and Doppler processor. 

Equation 42 has the same form that Equation 35 does. So after changing variables, 
setting kT s equal to zero and neglecting the first exponential term, the /th output of DFT 
can be written as 

Xi( x,f m ) =2 e^‘V(t-x,x)h*(-t)dt, (44) 

The above equation defines the cross-ambiguity function associated with the 
mathematical model of the stepped frequency radar processor as illustrated in Figure 10. 
The subscript l indicates that a distinct ambiguity function is defined for each Doppler 
filter. 

For the received pulse to be matched to the transmitted pulse, h(t) should be equal 
to j-i(-t). From Equations 37, 42, and 43, after simplifying, the ambiguity function for /th 
DFT can be expressed as 

min (N-l-l+pfi-l-IJrf-l) . i 

Xi(x,fd) = Z w(i) x e~ j2iu ^ x e-finV+'-p)^ x e flW<t-WKM)T r 

i=max (0^7-/) 


x(T s - iT r ) x x 


sm7t(/^my)(r^7V) 

7tfo-KA/)(2VyT r ) 


for 0 < y < Y r , 


(45-1) 
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min (N-I+p ft~I Jtf- 1) . l 

Xi(x,fd) = 2 w(i) x x x e /2 

j=max (0,/H-l-/) 


X f7\ + Cy- l'lTV') x e Mfd-W)( T s+ir-i)Tr) x )rr) 

^ s K ‘ > r) K(f d -KAf)(TM'f-l)Tr) ’ 


for 1 - j- < y < 1 , 


(45-2) 


elsewhere. 


(45-3) 


Taking the absolute value of the above equation, we obtain the ambiguity function 


= (TV -flV) x I lx 


mm (N-1 -/+/? jV-1 -/Jlf- 1) , 

2 vKO x e~ j2TU M x e jlK{fd ~ i2p ^ )Af)iTr I, 

f=max (0,/?-/) 

for 0 < y < ^, 


(46-1) 


(T i r- nr w I | y 

(J, + (Y i)ir)XI IX 


min (N-I+p fl-1 $1-1) . z 

2 vt<i) x e _ - /2ro w x eMfd-VmnWiTr | ; 

z=max (0^7+1-/) 

for 1 ~y <y< 1, 


(46-2) 


elsewhere. 


(46-3) 


The above equation becomes the cross-ambiguity function of the stepped 


frequency radar processor. The cross-ambiguity function can describe the auto-ambiguity 
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function of the stepped frequency radar by setting the Doppler filter index / to zero and 
assuming a uniform window, w(i)=l, setting / to zero, and setting the number of 
transmitted pulses, N, equal to the number of samples processed, M. Under these 
assumptions and after simplifying, the ambiguity function becomes 


\Xo(tJ d )\=(T s -ir r )x\ 


sMi{f<r-p£f){T s -ir r ) 

W<r-pAf)(T s -*(Tr) 


ixl 


sin Tt(fd-{2p^i)^)T r 


forO<y<£, (47-1) 

(T | (y i YT \ V I S ™^r~(^l)4/)(7>K , H)7V) I I sin7i(/j-(2jH^l)4/)(A f —ljH-ll)7V | 
t i U )*r) rifd-(p+VmTs+(r-l)Tr) 81 X 131 (^( 2 ^ 1 ) 4 /)^ ’ 



for l-£<y<l, 

(47-2) 

0 , 

elsewhere. 

(47-3) 


The above equation is exactly the same as the expression of the auto-ambiguity 
function in Equation 24. If we let the frequency step, A/, be zero, it is easy to prove that 
Equation 38 becomes the ambiguity function of a pulse train as expressed in Appendix A 
Equation. A. 13. 

The ambiguity diagram of the Oth, 1=0, Doppler filter of the stepped frequency 
radar with a Hamming window in the Doppler processor is illustrated in Figure 11. Figures 
12 and 13 are plots of cut along the frequency axis and delay axis respectively. From the 
comparison of Figures 12 and 13 with Figures 7 and 8 which are the cut plots of the 
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frequency (prf) 


-10 -4 


delay (pri) 


Figure 11. Cross-ambiguity Diagram of the Stepped Frequency Radar, Pulse Number 

N=4, Number of Pulses Processed M=4, Duty-cycle=0.2, A/ =2PRF, 1=0, /=0, 
and with Hamming Window 
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Figure 12. Plot of Cut along Frequency Axis for x = 0 of the Cross-ambiguity Function 
for the Stepped Frequency Radar, Pulse Number N=4, Number of Pulses 
Processed M=4, Duty cycle=0.2, A/=2PRF, 1=0, 1=0, and with Hamming 
Window 
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Figure 13. Plot of Cut along Time Axis for f=0 of the Cross-ambiguity Function for the 
Stepped Frequency Radar, Pulse Number N=4, Number of Pulses Processed 
M=4, Duty-cycle=0.2 ,A f =2PRF, 1=0, and 1=0, and with Hamming Window 
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auto-ambiguity diagram for the stepped frequency radar without window function, 
one can find that window functions which are often used in radar processors will decrease 
the amplitude of each spike and spread the width of it in return for reducing the sidelobes 
of each spike. The decrease in amplitude will reduce the detection capability and the 
spread of the width of the spike will cause poor accuracy and resolution. On the other 
hand, reducing the sidelobe of the spike will enhance the detection of small targets close to 
the large target. Figure 14 shows the bank of Doppler filters within a DFT processor . 
This figure indicates that the ambiguity function for the /th Doppler filter is the shifted 
version, along the frequency axis, of the ambiguity function for the Oth Doppler filter. The 
effect of decay for different Doppler filters is caused by the window function. 

The improvement in clutter suppression capability that is achieved using pulse 
Doppler processing, including the use of rejecting I pulses for initialization, is 
demonstrated in Figures 15 thru 28 [Ref. 7]. The figures were drawn using the MATLAB 
computer program which calculates the ambiguity function defined by Equation 47. 
Figures 15 thru 22 show the clutter suppression capability for the pulse radar, and Figures 
23 thru 28 show the clutter suppression capability for the stepped frequency radar. In all 
these plots, the sample step of time delay axis is 4 PRI. Therefore only envelopes of the 
ambiguity diagram are shown in the time delay axis. In the frequency shift axis, there are 
many spikes occurring at the position of multiples of PRF, but only a single spike is shown 
in the plots. The Hamming window is used in all of these plots. In Figures 15, 17, 19, 21, 
23, 25, and 28, the number of transmitted pulses and processed pulses are equal. 
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frequency (prf) 

Figure 14. Plot of Cut along Frequency Axis for t = 0 of the Cross-ambiguity Function 
for the Stepped Frequency Radar, Pulse Number N = 4 , Number of Pulses 
Processed M=4, Duty cycle= 0 . 2 , A/>2PRF, 1=0, /=[0, 1 , 2 , 3 ], and with 
Hamming Window 



N=M=100, and the number of initial pulses rejected, I, is set to zero. In Figures 16, 18, 
20, 22, 24, 26, and 28, the number of transmitted pulses, N, is increased to 150, and I is 
set to 50. From Figures 15 thru 18, in the case of the pulse radar, the use of initial pulse 
rejection leads to an ambiguity function which, along the delay axis, is a cross correlation 
of envelope functions of unequal length. This phenomenon is not clear for the stepped 
frequency radar due to the different pulse frequencies as shown in Figures 23 and 24. 
Comparing Figures 19, 21, 25, and 27 with Figures 20, 22, 26, and 28, respectively, the 
improvement in clutter suppression capability when utilizing initial pulse rejection is clearly 
demonstrated. Comparison of Figures 19 and 25 with Figures 20 and 26 shows that the 
clutter in the 0 to 50 PRI range will be rejected by the use of an initial pulse rejection for a 
0.1 PRF Doppler mismatch. Figures 21 and 27 show that, without initial pulse rejection, 
the capability of Doppler processing to suppress clutter is severely degraded when placing 
30 PRIs of delay (from the matched delay) when compared to the cases shown in Figures 
22 and 28. [Ref. 6] 
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Figure 15. Cross-ambiguity Diagram of the Pulse Radar, Pulse Number N=100, 
Number of Pulses Processed M=100, Duty-cycle=0.2, A/ =0,1=0, 1=0, 
and with Hamming Window 
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Figure 16. Cross-ambiguity Diagram of the Pulse Radar, Pulse Number N=150, 

Number of Pulses Processed M=100, Duty-cycle=0.2, A/ =0,1=50, 1=0, 
and with Hamming Window. 



Figure 17. Plot of Cut along Time Axis for f=0 of the Cross-ambiguity Function for the 
Pulse Radar, Pulse Number N=100, Number of Pulses Processed M=100, 
Duty cycle=0.2, A f= 0,1=0, /=0, and with Hamming Window. 



Pulse Radar, Pulse Number N=150, Number of Pulses Processed M=100, 
Duty-cycle=0.2, A/ =0,1=50, /=0, and with Hamming Window 
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Figure 19. Plot of Cut Parallel to Time Axis for f=0.02PRF of the Cross-ambiguity 


Function for the Pulse Radar, Pulse Number N=100, Number of Pulses 
Processed M=100, Duty-cycle=0.2, Af= 0,1=0, /=0, and with Hamming 
Window 



Figure 20. Plot of Cut Parallel to Time Axis for f=0.02PRF of the Cross-ambiguity 


Function for the Pulse Radar, Pulse Number N=150, Number of Pulses 
Processed M=100, Duty-cycle=0.2, A/=0,1=50, 1=0, and with Hamming 
Window. 
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Figure 21. Plot of Cut Parallel to Frequency Axis for t =30PRI of the Cross-ambiguity 
Function for thePulse Radar, Pulse Number N=100, Number of Pulses 
Processed M=100, Duty-cycle=0.2, A/= 0,1=0, 7=0, and with Hamming 
Window 



Figure 22. Plot of Cut Parallel to Frequency Axis for x =30PRI of the Cross-ambiguity 
Function for the PulseRadar, Pulse Number N=150, Number of Pulses 
Processed M=100, Duty-cycle=0.2, A f=0, 1=50, 1=0, and with Hamming 
Window. 
















Figure 23. Cross-ambiguity Diagram of the Stepped Frequency Radar, Pulse Number 
N=100, Number of Pulses Processed M=100, Duty-cycle=0.2 , Af =0.1PRF, 
1=0, /=0, and with Hamming Window 
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Figure 24. Cross-ambiguity Diagram of the Stepped Frequency Radar, Pulse Number 
N=150, Number of Pulses Processed M=100, Duty-cycle=0.2 , A/ =0.1PRF, 
1=50, 1=0, and with Hamming Window 































Figure 25. Plot of Cut Parallel to Time Axis for f=0.02PRF of the Cross-ambiguity 

Function for the Stepped Frequency Radar, Pulse Number N=100, Number of 
Pulses Processed M=100, Duty-cycle=0.2, A/’=0.1PRF, 1=0, 1=0, and with 
Hamming Window 



Figure 26. Plot of Cut Parallel to Time Axis for f=0.02PRF of the Cross-ambiguity 


Function for the Stepped Frequency Radar, Pulse Number N=150, Number of 
Pulses Processed M=100, Duty-cycle=0.2 ,A/’=0.1PRF, 1=50, 7=0, and with 
Hamming Window. 
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Figure 27. Plot of Cut Parallel to Frequency Axis for t=30PRI of the Cross-ambiguity 

Function for the Stepped Frequency Radar, Pulse Number N=100, Number of 
Pulses Processed M=100, Duty-cycle=0.2, A/'=0.1PRF, 1=0, 7=0, and with 
Hamming Window. 



Figure 28. Plot of Cut Parallel to Frequency Axis for x=30PRI of the Cross-ambiguity 

Function for the Stepped Frequency Radar, Pulse Number N=150, Number of 
Pulses Processed M=100, Duty-cycle=0.2, A/ =0.1PRF, 1=50, 7=0, and with 
Hamming Window. 
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IV. CONCLUSIONS 


This thesis investigates the ambiguity functions of the stepped frequency waveform 

and the stepped frequency radar system. The stepped frequency waveform achieves a high 

range resolution by coherently processing the returns from N pulses, each having a 

different earner frequency that changes by a fixed amount from pulse to pulse. The 

mathematical expression of the auto-ambiguity function for the stepped frequency 

waveform is derived from the definition. The auto-ambiguity function is defined as the 

correlation of the complex envelope of the transmitted signal and the Doppler frequency 

shifted version of the complex conjugate of the transmitted signal. The 3D plot of the 

auto-ambiguity function shown in Figure 5 is spiky like that of a constant frequency pulse 

train ( shown in Figure A.5) and is tilted at an angle. Along the time delay axis, the 

null-to-null width of each spike is 2/N Af as compared to 2/x for the constant frequency 

pulse train, thereby making it possible to decrease the width by increasing NA/ without 

increasing its instantaneous bandwidth. However, for the constant frequency pulse train, 

high resolution can be achieved by decreasing the pulse width and, thereby, increasing the 

instantaneous bandwidth. On the other hand, the null-to-null spike width of the stepped 
frequency waveform along the frequency axis is equal to 2/(NxPRI) which is the same as 

for the constant frequency pulse train. Thus frequency resolution is not improved by the 

stepped frequency waveform. 
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The mathematical expression of the cross-ambiguity function for the high PRF 
stepped frequency radar system including its receiver and signal processor is also derived. 
This mathematical expression is obtained from the output of the Doppler processor and a 
distinct ambiguity function is defined for each Doppler processing filter. If weighting is 
ignored, the cross-ambiguity function of the Oth Doppler filter reduces to the 
auto-ambiguity function of the stepped frequency waveform derived from the traditional 
definition. Ambiguity functions of other Doppler filters are shifted versions of the 
auto-ambiguity function along the frequency axis. The mathematical expression has been 
verified by comparing the ambiguity diagram obtained by simulation which performs the 
correlation method. 

The improvement of the clutter suppression capability for the stepped frequency 
radar system by rejecting initial pulses is also proven. The strong clutter returns from close 
range can be suppressed by rejecting some initial pulses, thereby reducing the interference 
from the region. Four types of basic waveforms are also discussed in the appendix for the 
purpose of comparison. These are single pulse, periodic constant frequency pulse train, 
linear frequency modulated pulse, and discrete frequency modulated pulse. The 
mathematical expressions of ambiguity functions for these waveforms are derived and 
have been also verified by simulation. 
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APPENDIX A. EXAMPLES OF THE AMBIGUITY FUNCTION FOR THE 
BASIC WAVEFORMS 

In this section the ambiguity functions of commonly used radar signals such as 
single pulse, constant frequency pulse train, linear frequency modulated pulse, and discrete 
frequency modulated pulse are discussed. The mathematical expression of the ambiguity 
function for these waveforms are derived and verified by comparing with the ambiguity 
diagram obtained from the simulation. 

1. Single Pulse of the Sine Wave 

The single pulse of the sine wave can be defined as 

St(t) = u(t)e^ , (A.l) 

where u(t) is the complex envelope of the signal, it is defined as 


u(t) = 1, if 0< t<T s , 


= 0, elsewhere, and 

T s is the pulse width. 


(A.2) 

The ambiguity function of the single pulse can be expressed as 
X(x,fd) = J u(t—x)u*(t)e j2 ^ dt dt 

—OO 


T 

=| e j2r ^ dt dt 

for 0 < x < T s , 

(A.3-1) 

X(T M - (J, + X) x X , 

for —T s < x < 0, 

(A.3-2) 

X(x,f d ) = 0, 

elsewhere. 

(A.3-3) 
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The absolute value of the ambiguity function of the single pulse can be 


expressed as 


lw*)l=(n-T)x|2gg2 


, for 0 < x < T s , 


(A4-1) 


iw*)i=<r.+T)xi**iga2i, f„ r - r ,< x <o, (A4-2) 

=0, elsewhere. (A4-3) 

The ambiguity function along the Doppler frequency shift axis, x=0, 

becomes 


\X(0,f d )\=T s x\^^\ 


(A. 5) 

The ambiguity function along the time delay axis,/ d =0, becomes 


lx(x,o)l=r,-x , 

for 0 <x<T s , 

(A.6-1) 

= T s + t , 

for -T s <x<0, 

(A.6-2) 

=0, 

elsewhere. 

(A.6-3) 


The ambiguity diagram, IX(x/ d )l 2 , of the single pulse is shown in Figure 
A. 1, the contour plot is shown in Figure A.2, and plots of cut along the frequency axis for 
x=0 and time delay axis for / d =0 are shown in Figure A.3 and Figure A.4. The plot for 
zero time delay is the square of the spectrum of a rectangular pulse which is the square of 
the absolute value of the SINC function. The plot for zero frequency shift represents the 
magnitude square of the auto-correlation function of a rectangular pulse. 
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Figure A. 1. Ambiguity Diagram of A Single Pulse, Pulse Width=l. 
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Figure A.2. Contour Plot of the Ambiguity Diagram of A Single Pulse, Pulse Width=l. 
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frequency (1/pu1se - width) 


Figure A.3. Plot of Cut along Frequency Axis for t=0 of the Ambiguity Function for 
A Single Pulse, Pulse Width=l. 



Figure A. 4. Plot of Cut along Time Axis for f=0 of the Ambiguity Function for 
A Single Pulse, Pulse Width=l. 
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2. Constant Frequency Pulse Train 


Constant frequency pulse train of N pulses can be defined as 


S t (t) = Z ]L(t-nT r 

n=0 


(A.7) 


where p(t) , the complex envelope of the single pulse of transmitted signal, is defined as 


where 


\x(t-nT r ) = 1, 

0 , 


if nT r ^t< nT r + T s , 
elsewhere, 


(A. 8) 


N : the number of pulses, 

77 : pulse repetition interval (PRI), 
fo : carrier frequency. 

The ambiguity function of periodic pulse train can be derived as 


X(x,f d ) = J E 1 ]±(t-nT r -x) E 1 ( x*(t-mT r )eJ 2 ^dt 

—~n=0 m=0 


(A.9) 


After changing variable, we obtain the ambiguity function as 


V- IN-1 00 

X(zJ d ) = E E eW* nTr J p(r)p* n)T r - x)e^“ l dt 

n=0m=0 —OO 


(A. 10) 


Now we define the round trip time in terms of integer and fraction parts 
x = (p + y)T r (A. 11) 


where p is integer ,y is the fraction, and 0 < y < 1. 
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After substituting Equation A. 11 in Equation A. 10 and simplifying , we 


obtain the expression for the ambiguity function as 


X(x,f d ) = C T s -iT r ) x x sm J^ Ts ~l T ; ) x^ {N Y~ 1+P) e^ nTr , 

KfATs^Tr) n=max (0> p) 


for 0<y<y, 


(A.12-1) 


( T s + (y-1 )T r ) x ^(A-Kr-Drr) x ^„ Tr , 

r) «=max (0^H-1) 


for 1 <y< 1, 


(A. 12-2) 


elsewhere. 


(A. 12-3) 


We can now obtain the ambiguity function as the following equation after 
simplifying the summation term and taking the absolute value. 


lv/_ / \ l _/T a rr \ ^\^d(Js^Tr)\ _ I smvf d (N-\p\)T, | 

\X(Z,f d )\ - (Ts-iTr) x I ■^Jr'^Tr) 1 X 1 Sin Wr 1 ’ 


for 0<y<|f, 


(A13-1) 


lv/_ ,rr , i \t \ „, | Sin^(^+(7-1)1,) | w | smvf d (N-\p+W)T r | 

IZ(T,/ d )l - (Ts + (y~ Wr) X \-^ d(TMl _, )Tr) I X I--1 , 

for 1 - y < y < 1, 


(A 13-2) 


\X(x,f d )\=0, 


elsewhere. 


(A13-3) 
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The ambiguity function along the frequency axis, i.e., T = 0, can be 


expressed as 


\X(0,f d )\=T s x\ 


si &WdTs) 
TfdTs 


sm(nNf d T r ) 

sm(nf d Tr) 


(A. 14) 


as 


The ambiguity function along the time delay axis, i.e., fd = 0, is expressed 

te,o)Uov-lpl)(r,-YT r ), 

for 0 <y< (A.15-1) 


=(jv-lp+il)((Y-i)T r +r f ), 

for l-|f <y< 1, (A. 15-2) 

= 0 , elsewhere. (A. 15-3) 

The ambiguity diagram of the pulse radar according to Equation A. 13 is shown in 

Figure A.5 and the contour of the ambiguity diagram is shown in Figure A.6. Plots of cut 

along frequency axis for x=0 and cut along time delay axis for/ d =0 are shown in Figure 

A.7 and A.8. From Equation A. 14 and Figure A.7, spikes appear at frequencies which are 

multiples of PRF. The first null for each spike is l/NT r away from the center of the spike. 

Therefore, each spike has the null to null width of 2/AT) in the frequency axis. Amplitudes 

of those spikes varies as a SINC function. Zero amplitude will appear at 1/7). From 

Equation A. 15 and Figure A.8, it is seen that there are 2N-1 spikes along the time delay 

axis each of which has the width of 27). The plot of cut of the ambiguity function at zero 

frequency shift represents the magnitude square of the auto-correlation function of a pulse 

train. 
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Figure A.5. Ambiguity Diagram of the Pulse Radar, Number of Pulses N=4, 
Duty-cycle=0.2. 
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Figure A.6. Contour Plot of the Ambiguity Diagram for the Pulse Radar, Number of 
Pulses N=4, Duty-cycle=0.2. 
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Figure A.7. Plot of Cut along Frequency Axis for t=0 of the Ambiguity Function for the 
Pulse Radar, Number of Pulses N=4, Duty-cycle=0.2. 



Figure A. 8. Plot of Cut along Time Axis for f=0 of the Ambiguity Function for the Pulse 
Radar, Number of Pulses N=4, Duty-cycle=0.2. 
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3. Linear Frequency Modulated Pulse 

Linear frequency modulated pulse can be mathematically represented as 




= [uCOe^y 2 ^* t 


(A.16) 


where u(t) was defined in Equation A.2 and |i is the rate of frequency change. The 
ambiguity function of the linear frequency modulated pulse can be written as 


X(zJ d ) =j u(t— x ) e wM 2 u * (Oe-^e^dt 




(A. 17) 


Equation A. 17 is exactly the same as Equation A.3, except that/ d is 


replaced by/ d -|ux. The ambiguity function becomes 


\x(x,u) I =(r, - x) x I 1 , forost<r„ 


(A.18-1) 


for-r^TcO. (A. 18-2) 


elsewhere. 


(A. 18-3) 
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The ambiguity diagram is shown in Figure A.9 and contour plot is shown 
in Figure A. 10. From Figure A.2 and Figure A. 10, it can be seen that the ambiguity 
function of the LFM pulse is a rotated version of the ambiguity function of a single pulse. 
The angle of rotation depends upon the rate of frequency change |i as shown in Figure 
A.11. 

4. Discrete Frequency modulated Pulse 

The discrete frequency modulated pulse consists of N subpulses, each at 
different discrete carrier frequency. Such a pulse can be expressed as 

S t (t) =2 u(t - n T)e J1 ^ +n ^ t , (A. 19) 

n=0 

where N is the number of the discrete frequencies in the pulse, NT is the pulse width, and 
u(t) is defined as 

u(t-nT) = 1, if nT<t<(n+l)T, 

= 0, elsewhere. (A.20) 

The envelope of the transmitted signal can be written as 

U(t) = X u(t - nT)e j2m£ ^. (A.21) 
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Figure A.9. Ambiguity Diagram of the LFM Pulse, Pulse Width=l. 
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Figure A. 10. Contour Plot of the Ambiguity Diagram of the LFM Pulse, Pulse Width=l. 
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Figure A. 11 .The Relation of the Ambiguity Functions for the Single Pulse and LFM Pulse. 
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According to the definition of the ambiguity function, we obtain the 


ambiguity function of the discrete frequency modulated pulse as 


X(x,f d )=j U(t-'L)U*{t)eW“ t d\. (A.22) 


Now we define the round trip time in terms of integer and fraction parts, 

x = (p+y)7\ (A.23) 

where p is an integer and y is a fraction where 0 < y < 1. We can obtain the ambiguity 

function of the discrete frequency modulated pulse after simplification. The expression of 
the ambiguity function has the form as 


X(x,fd) = {(1 - y)T x e Md-p&f)(i+y)T Xe j2np&fz x e jK(fd-(2p+yW)(N-i+p)T x 


sinirQ^-pA/Kl^r sm7c(/^(2p+7)A/)(A'-l;>l)7' ' 
siMt(f d -(2pvf)&f)T -* 

{y J x e M<Hp+Y)bf)1T-xpn(p+V)6fi x e M<r-(2p+y+l)W)(N+p)T x 


^nTUfj-Xp+lWpfT 

K(f d -(p+l)Af)yT 


sin Ti(f d ~(2p+y+l)&f)(N-\p+l\)T 
sin it(f cr -(2p+yf l)A/)r 


}• 


(A.24) 


The ambiguity diagram is shown in Figure A. 12 and contour plot is shown 
in Figure A. 13. The discrete frequency modulated pulse is composed of the linear 
frequency modulated signal and the pulse train with 100% duty-cycle. Therefore, the 
ambiguity diagram of it has been tilted and has many discrete spikes. 
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Figure A. 12. Ambiguity Di, 
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Figure A. 13. Contour Plot of the Ambiguity Diagram of the Discrete FM Pulse, 
Pulse Width=l. 





APPENDIX B. AMBIGUITY FUNCTION PROGRAM CODES 


%filename: progl.m 

% Ambiguity function of the single pulse by correlation method 
%update: 27/sep/94 
%by: Huang, Jen-chih 

% 

clear; 
tau=l; 
nx=50; 
ny=101; 

t=[0:nx-l]*tau/nx; 
dt=t(2)-t(l); 

z=0; 

fd=linspace(-4/tau,4/tau,ny) 
ul=ones(l,nx); 
for m=l:ny 

u2=u 1. *exp(j *2*pi*fd(m) *t); 
c=xcorr(u2,ul).*dt; 
z=[z;(abs(c)). A 2]; 
end 

t=[fliplr(-t),t(2:nx)]; 
figure(l) 
mesh(t,fd,z) 
title(’Ambiguity diagram of single pulse') 
xlabel('delay (pulse-width)') 
ylabel('frequency (1/pulse-width)') 
figure(2) 

plot(t,z((ny+l)/2,:)) %plot of cut along frequency axis 

xlabel('delay (pulse-width)') 
ylabel('frequency (1/pulse-width)') 
figure(3) 

plot(fd,z(:,nx)) %plot of cut along delay axis 

xlabelfdelay (pulse-width)') 
ylabel('frequency (1/pulse-width)') 


% clear all the variable 
%pulse width =1 
%# of step in delay axis 
%# of step in frequency axis 
%time axis 
%time step 

%frequency axis 

%envelope of the transmitted signal 

%ul multiply by doppler freq. shift 
% correlation of ul & u2 
%matrix of the ambiguity function 

%reconstruct the time axis 
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^filename: prog2.m 

% Ambiguity function of the linear frequency modulation (LFM) pulse 
%by correlation method 
% update: 15/aug/94 
%by: Huang, Jen-chih 

% 

% 

clear; 
tau=l; 
mu=4; 
nx=50; 
ny=101; 

t=[0:nx-l]*tau/nx; 
dt=t(2)-t(l); 
z=0; 

fd=linspace(-4/tau,4/tau,ny); 
for m=l:ny 

ul=exp(j*pi*mu.*t*t); 
u2=conj(ul).*exp(j*2*pi*fd(m)*t); 
c=conv (u2,fliplr(u 1)). *dt; 
z=[z;(abs(c)). A 2]; 
end 
% 

t=[fliplr(-t),t(2:nx)]; %reconstruct the time axis 

figure(l) 

mesh(t,fd,z) %mesh plot of ambiguity function 

title('Ambiguity function of the single pulse') 

xlabel('delay (pulse-width)') 

ylabel('frequency (1/pulse-width)') 

figure(2) 

contour(t,fd,z,40) %contour plot of the ambiguity function 

title('Contour plot of the signal pulse') 
xlabel('delay (pulse-width)') 
ylabel('frequency (1/pulse-width)') 

% 


%clear all the variables 

%pulse width=l 

%rate of frequency change 

%# of step in delay axis 

%# of step in frequency axis 

%time axis 

%time step 

%frequency axis 

%envelope of the transmitted signal 
%ul multiply Doppler frequency shift 
%correlation of ul & u2 

%matrix of the mag. square of the ambiguity function 
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% filename: prog3.m 

% Ambiguity function of the stepped frequency radar by the correlation method 
%update: 15/mar/94 
%by: Huang, Jen-chih 
% 


clear; 


tau=0.2; 

t0=0.8; 

df=0; 

% 

% 

fs=40; 

N=4; 

T=N*(tau+tO); 

sn=fs*tau; 

tn=fs*tO; 

Tn=sn+tn; 

P=[ones(l,sn),zeros(l,tn)]; 

t=linspace(0,T,N*Tn); 

nx=101; 

fd=linspace(-10,10,nx); 
dt=t(2)-t(l); 

s=0; 


%pulse width 

%PRI-pulse-width; PRI=1 
%step frequency 
% if df=0 then pulse radar 
%else stepped frequency radar 
% sample rate (Hz) 

%# of pulses 

%total time period for N pulses 
%# of samples in pulse width 
%# of samples between pulses 
% of samples in one PRI 
% envelope for one PRI 
%time axis for N pulses 
%# of samples in frequency axis 
%frequency axis 
%time step 


m=[l,2,3,4]; %step frequency sequence 

for n=l:N 

s=[s,exp(j*2*pi*(m(n)-l)*df*t(Tn*(n-l)+l:Tn*n)).*P]; %transmitted signal 
end 


z=[|; 

for L=l:mx 
f=fd(L); 

s 1 =conj (s).*exp(j *2*pi*f*t); 
s2=s; 

c=conv (s 1 ,fliplr(s2)).*dt; 

z=[z;(abs(c)). A 2]; 

end 

[l,r]=size(z); 
t=[fliplr(-t),t(2:N*Tn)]; 
figure(l) 
mesh(t,fd,z) 
xlabel('delay (PRI)') 
ylabel('ffequency (PRF)') 


%transmitted signal multiply by doppler shift 
%the transmitted signal 
%correlation of si & s2 
%matrix of the ambiguity function 


%resconstruct time axis 
%mesh plot of the ambiguity function 
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figure(2) 

plot(t,z((nx+l)/2,:)) %plot of cut along time axis 

xlabel('delay (PRI)’) 

grid 

figure(3) 

plot(fd',z(:,(r+l)/2)) %plot of cut along frequency axis 

xlabel('frequency (PRF)') 

grid 

% 
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*5fc ***** ** * * * *** ** * *** * ****** * ******* **** ** 5jc»fc^c5fc5fc5fc5fc5fc *** 5jc2f= * ********** ** *** * 


%filename: prog4.m 

%Ambiguity function of the discrete frequency pulse 
%update: 18/max/94 
%by: Huang, Jen-chih 

% 

% 

clear %clear all the variables 


tau=l; 

df=l; 

N=4; 

sn=30; 


ny=101; 

t=[0;N*sn-l]*tau/(N*sn); 

dt=t(2)-t(l); 

s=D; 

m=[l,2,3,4]; 
for n=l:N 


%pulse width 
%frequency step 
%# of discrete frequencies 
% sample rate in delay axis 
%# of sample in frequency axis 
% time delay axis 
%time step 

%discrete frequency sequence 


s=[s,exp(j*2*pi*(m(n)-l)*df*t(sn*(n-l)+l:sn*n))]; %envelope of the transmitted signal 
end 


z=[]; 

fd=linspace(-4/tau,4/tau,ny); %ffequency axis 

for 1=1 :ny 

sl=conj(s).*exp(j*2*pi*fd(l)*t); %conjugate of s times Doppler shift 

s2=s; 


c=conv(sl,fliplr(s2)).*dt; %correlation of si & s2 

z=[z;(abs(c)). A 2]; %matrix of the ambiguity function 

end 


t=[fliplr(-t),t(2:N*sn)]; %reconstruct the time delay axis 

figure(l) 

mesh(t,fd,z) %mesh plot of the ambiguity function 

title('Ambiguity function of the discrete frequency pulse') 

xlabel(' time (pulse width)') 

ylabel('frequency (1/pulse-width)') 

figure(2) 

contour(t,fd,z,40) %contour plot of the ambiguity function 

title('contour plot of the ambiguity function of the discrete frequency pulse') 
xlabel(' time (pulse width)') 
ylabel('frequency (1/pulse-width) 1 ) 

% 


* * * * * * * * * * * * ** ** * * * * * *** *** * ** * * * * ** * *** * ** * * * * * * *** ** ** **** ** * * * * * * * ** * 
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%filename: progal.m 

%ambiguity function of the single pulse by equation 
%update: 20/aug/94 
%by: Huang, Jen-chih 

% 

% 

clear; 

Ts=l; 

nx=101; 

tau=linspace(-Ts,T s,nx); 
ny=51; 

fd=linspace(-4/Ts,4/Ts,ny); 

z=D; 

for m=l :ny 
f=fd(m); 

x=[]; 

for n=l:nx 
t=tau(n); 
if t>=0 

if f*(Ts-t)==0 
x=[x,(Ts-t)]; 
else 

x=[x,(Ts-t). * (sin(pi*f*(T s-t))) ./(pi*f* (T s-t))]; 
end 
else 

if f*(Ts+t)==0 
x=[x,(Ts+t)]; 
else 

x=[x, (T s+t). *(sin(pi*f*(T s+t)))./(pi*f*(T s+t))]; 
end 
end 
end 

z=[z;x]; 

end 

z=(abs(z)). A 2; 
mesh(tau,fd,z) 
xlabel('delay (pulse-width)') 
ylabel(’frequency (1/pulse-width)') 
contour (tau,fd, z,40) 

sX- si# sl> \X- \L# sL* si# si# s’# sL* *X 'X - 4* *1# si# si# si# si# si# si# si# -X- si# *1# si# X* #1# sZ# s’# sL# sl> sl> si# si# si# si# si# si# si# si# si# si# si# si# si# si# si# si# si# si# si# si# si# si# si# si# si# si# si# si# 
•T* ■'T* "T* v ^ ‘*T* ■'T* •T* v ■'T* «T* "T* •T* T 1 *T“ v v *7* •T* *1* v 'T* "T* ^T* -T* 'T* «T* *£* V #r* ^ <p -p /f* ip <p ip ip #p #p ip #p <p ip 
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%filename: proga2.m 

% Ambiguity function of the linear frequency modulation pulse by equation 

% 

clear; 

Ts=l; 

mu=4; 

nx=101; 

tau=linspace(-T s,T s,nx); 
ny=51; 

fd=linspace(-4/T s ,4/T s,ny); 

z=D; 

for m=l:ny 
x=[]; 

for n=l:nx 
t=tau(n); 
f=fd(m)-mu*t; 
if t>=0 

if f*(Ts-t)==0 
x=[x,(Ts-t)]; 
else 

x=[x,(T s-t). * (sin(pi*f* (T s-t)))./(pi*f* (T s-t))]; 
end 
else 

if f*(Ts+t)=0 
x=[x,(Ts+t)]; 
else 

x=[x, (T s+t). * (sin(pi*f* (T s+t)))./(pi*f* (T s+t))]; 
end 
end 
end 

z=[z;x]; 

end 

z=(abs(z)). A 2; 

mesh(tau,fd,z) 

xlabel('delay (pulse-width)') 

ylabel('frequency (1/pulse-width)') 

contour(tau,fd,z,40) 

% 
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%filename: proga3.m 

% Ambiguity function of the pulse radar by equation 
% update: 20/aug/94 
%by: Huang,Jen-chih 
clear; 

Ts=0.2; 

Tr=l; 

d=Ts/Tr; 

N=4; 

nx=161; 

tau=linspace(-N*Tr,N*Tr,nx); 

ny=61; 

fd=linspace(-10,10,ny); 

p=floor(tau/Tr); 

r=(tau/Tr)-p; 

z=D; 

for m=l:ny 
x=[]; 
f=fd(m); 
for n=l:nx 
t=tau(n); 
pp=p(n); 
rr=r(n); 

if rr>=0 & rr<d 
if f==0 

x=[x,(Ts-rr*Tr)*(N-abs(pp))]; 

else 

x 1 =(sin(pi*f* (T s-rr*Tr)))./(pi*f* (T s-rr*Tr)); 
if rem(f*Tr,l)==0 
x2=N-abs(pp); 
else 

x2^sin(pi*f*(N-abs(pp))*Tr))V(sin(pi*f*Tr)); 

end 

x=[x,(Ts-rr*Tr).*xl*x2]; 
end 
else 
if f=0 

x=[x,(Ts+(rr-l)*Tr)*(N-abs(pp+l))]; 

else 

xl=( sin (pi*f* (T s+(rr-l)*Tr)))./(pi*f* (T s+(rr-1 )*Tr)); 
if rem(f*Tr,l)==0 
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x2=N-abs(pp+l); 

else 

x2=(sin(pi*f*(N-abs(pp+l))*Tr))./(sin(pi*f*Tr)); 

end 

x=[x,(Ts+(rr-l)*Tr).*xl*x2]; 

end 

else 

x=[x,0]; 

end 

end 

z=[z;x]; 

end 

z=(abs(z)). A 2; 

figure(l) 

mesh(tau,fd,z) 

xlabel('delay (pulse-width)') 

ylabelCfirequency (1/pulse-width)') 

figure(2) “ 

contour(tau,fd,z,40) 

% 

************************************************************************ 
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%filename: proga4.m 

% Ambiguity function of the stepped frequency radar by equation 
% update: 20/aug/94 
%by: Huang, Jen-chih 
clear; 

Ts=0.2; 

Tr=l; 

df=2; 

d=Ts/Tr; 

N=4; 

nx=161; 

tau=linspace(-N*Tr,N*Tr,nx); 

ny=61; 

fd=linspace(-10,10,ny); 

p=floor(tau/Tr); 

r=(tau/Tr)-p; 

z=Q; 

for m=l:ny 

x=[]; 

f=fd(m); 
for n=l:nx 
t=tau(n); 
pp=p(n); 
rr=r(n); 
fl=f-pp*df; 
f2=f-(2*pp+rr)*df; 
f3=f-(pp+l)*df; 
f4=f-(2*pp+rr+1) *df; 
if rr>=0 & iT<d 
if f 1 ==0 
xl=l; 
else 

xl=(sin(pi*fl*(Ts-rr*Tr)))./(pi*fl*(Ts-rr*Tr)); 

end 

if rem(f2*Tr,l)==0 
x2=N-abs(pp); 
else 

x2=(sin(pi*f2*(N-abs(pp))*Tr))./(sin(pi*f2*Tr)); 

end 

x=[x,(Ts-rr*Tr).*xl*x2]; 
elseifrr<l &rr>=(l-d) 
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if £3=0 
xl=l; 
else 

x 1=(sin(pi*f3*(T s+(rr-1 )*Tr)))./(pi*f3*(T s+(rr-1 )*Tr)); 
end 

if rem(f4*Tr,l)==0 
x2=N-abs(pp+l); 
else 

x2=(sin(pi*f4*(N-abs(pp+l))*Tr))./(sin(pi*f4*Tr)); 

end 

x=[x, (Ts+(rr-1) *Tr). *x 1 *x2]; 

else 

x=[x,0]; 

end 

end 

z=[z;x]; 

end 

z=(abs(z)). A 2; 

figure(l) 

mesh(tau,fd,z) 

xlabel('delay (pulse-width)') 

ylabel('frequency (1/pulse-width)') 

figure(2) 

contour(tau,fd,z,40) 

% 
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^filename: proga5.m 

% Ambiguity function of the discrete frequency modulation pulse by equation 
% update: 20/aug/94 
%by: Huang, Jen-chih 
clear; 

N=4; 

T=l/N; 

df=2; 

d=Ts/Tr; 

nx=121; 

tau=linspace(-N*Tr,N*Tr,nx); 

ny=61; 

fd=linspace(-4,4,ny); 

p=floor(tau/Tr); 

r=(tau/Tr)-p; 

z=D; 

form=l:ny 
x=[]; 
f=fd(m); 
forn=l:nx 
t=tau(n); 
pp=p(n); 
rr=r(n); 
fl=f-pp*df; 
f2=f-(2*pp+rr)*df; 
f3=f-(pp+l)*df; 
f4=f-(2*pp+rr+l)*df; 
if fl*(Ts-rr*Tr)==0 
xl=l; 
else 

xl=(sin(pi*fl*(Ts-rr*Tr)))./(pi*fl*(Ts-rr*Tr)); 

end 

if rem(f2*Tr,l)==0 
x2=N-abs(pp); 
else 

x2Ksin(pi*f2*(N-abs(pp))*Tr))V(sin(pi*f2*Tr)); 

end 

yl=(Ts-rT*Tr).*exp(j*pi*fl*(Ts+rr*Tr)).*xl.*exp(j*2*pi*pp*df*t) 

. *exp (j *pi *f2*(N-1 +pp) *Tr). *x2; 
if f3*(Ts+(rr-l)*Tr==0 
x3=l; 
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else 

x3=(sin(pi*f3 * (T s+(rr-1) *Tr)))./(pi*f3*(T s+(rr-1 )*Tr)); 
end 

if rem(f4*Tr,l)==0 
x4=N-abs(pp+l); 
else 

x4=(sin(pi*f4*(N-abs(pp+l))*Tr))./(sin(pi*f4*Tr)); 

end 

y2=(Ts+(rr-l)*Tr).*exp0*pi*f3*(Ts+(rr-l)*Tr)).*x3.*exp(j 5f: 2*pi*(pp+l)*df i: t) 
. *exp (j *pi*f4*(N+pp) *Tr). *x4; 

x=[x,yl+y2]; 

end 

z=[z;x]; 

end 

z=(abs(z)). A 2; 

figure(l) 

mesh(tau,fd,z) 

xlabel('delay (pulse-width)') 

ylabel('frequency (1/pulse-width)') 

figure(2) 

contour(tau,fd,z,40) 

% 
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